Cluster form factor calculation in the ab initio no-core shell model 
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We derive expressions for cluster overlap integrals or channel cluster form factors for ab initio 
no-core shell model (NCSM) wave functions. These are used to obtain the spectroscopic factors and 
can serve as a starting point for the description of low-energy nuclear reactions. We consider the 
composite system and the target nucleus to be described in the Slater determinant (SD) harmonic 
oscillator (HO) basis while the projectile eigenstate to be expanded in the Jacobi coordinate HO 
basis. This is the most practical case. The spurious center of mass components present in the SD 
bases are removed exactly. The calculated cluster overlap integrals are translationally invariant. As 
an illustration, we present results of cluster form factor calculations for ( 5 He| 4 He+n), ( 5 He| 3 H+d), 
< 6 Li| 4 He+d), ( 6 Be| 3 He+ 3 He>, ( 7 Li| 4 He+ 3 H>, ( 7 Li| 6 Li+n), < 8 Be| 6 Li+d), ( 8 Be| 7 Li+p), ( 9 Li| 8 Li+n) 
and ( 13 C| 12 C+n), with all the nuclei described by multi-7if2 NCSM wave functions. 

PACS numbers: 21.60.Cs, 21.45.+V, 21.30.-x, 21.30.Fe 



I. INTRODUCTION 



There has been a significant progress in the ab initio approaches to the structure of light nuclei. Starting from the 
' realistic two- and three-nucleon interactions the methods like the Green's Function Monte Carlo (GFMC) [lj or the 
ab initio no-core shell model (NCSM) can predict the low-lying levels in p-shell nuclei. It is a challenging task to 
extend the ab initio methods to describe nuclear reactions. This is in particular true for low-energy reactions where 
detailed knowledge of nuclear structure is important. The first capture reaction calculations using the GFMC (or 
rather variational Monte Carlo) wave functions were performed 2] . Concerning the NCSM, in order to make the first 
steps in this direction one needs to understand the cluster structure of the eigenstates, i.e. to calculate the channel 
cluster form factors. Those can then, e.g., be integrated to obtain the spectroscopic factors. At the same time, starting 
from the channel cluster form factors one can attempt to set up an approach in the spirit of the resonating group 
method (RGM) to calculate radial wave functions describing the relative motion of the binary clusters and then 
q ■ obtain the cross sections. This paper addresses the issue of channel cluster form factor (or cluster overlap integral or 
^5 , reduced width amplitude for two-body decay) calculations in the NCSM. 

The principal foundation of the ab initio NCSM approach is the use of effective interactions appropriate for the 
large but finite basis spaces employed in the calculations. These effective interactions are derived from the underlying 
realistic inter-nucleon potentials through a unitary transformation in a way that guarantees convergence to the exact 
solution as the basis size increases. For the basis, one uses antisymmetrized A-nucleon harmonic-oscillator (HO) 
states that span the complete N max h£l space. A disadvantage of the HO basis is its unphysical asymptotic behavior, 
a problem that must be dealt with by using a large basis expansion and/or a renormalization. On the other hand, the 
nuclear system is translationally invariant and, in particular in the case of light nuclei, it is important to preserve this 
symmetry. The HO basis is the only basis that allows a switch from Jacobi coordinates to single particle Cartesian 
coordinates without violating the translational invariance. Consequently, one may choose the coordinates according 
to whatever is more efficient for the problem at hand. In practice, it turns out that the A = 3 system is the easiest 
solved in the Jacobi basis, the A = 4 system can be solved either way with the same efficiency when only two-body 
interaction is utilized, but the Jacobi basis is more efficient when the three-body interaction is included. For systems 
with A > 4, it is by far more advantageous to use the Cartesian coordinates and the Slater determinant (SD) basis 
and employ the powerful shell model codes like Antoine |]| that rely on the second quantization techniques. While 
the NCSM eigenenergies are independent on the choice of coordinates, the eigenfunctions obtained in the Cartesian 
coordinate SD basis include a Ohfl spurious center of mass (CM) component. 

Our goal is to calculate the channel cluster form factors regardless of the choice of coordinates. Obviously, the most 
desired case is the one corresponding to the most efficient choice, i.e., the projectile, that is the lighter nucleus of 
the binary system, consisting of a < 4 nucleons described by a Jacobi coordinate wave function, while the A-nucleon 
composite system and the A — a nucleoli target, that is the heavier nucleus of the binary system, described by wave 
functions expanded in the SD basis. To obtain the physical, translationally invariant cluster form factors we must 
remove completely the spurious CM components. 

The ways how to remove these c omp onents and obtain physical matrix elements of different operators were investi- 
gated in the past |rl PI lil 1^. ITol ITU Il2f . At the same time, the single-nucleon as well as cluster overlap integral and /or 
spectroscopic factor calculations were investigated in many papers, see e.g. Refs. [T3L [T^. f]~5L IT^. TT^ . Tl3 . IT^ . l20l . I2TL . 
In many cases, however, the basis space was limited to a single major HO shell. In the NCSM, the basis space spans 
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several major shells. In general, it is necessary to re- visit and adapt the techniques of the channel cluster form factor 
and spectroscopic factor calculation and the spurious center-of-mass motion removal to make them applicable for the 
NCSM. In an earlier investigation, we addressed the spurious center-of-mass motion problem for the density operator 
[23|. In this paper, we focus on the calculation of the channel cluster form factors. 

In Section^ we present the derivation and the algebraic formulas for calculating the channel cluster form factors 
from the NCSM wave functions for projectiles consisting of up to three nucleons. In Section lTTTl applications to several 
light nuclei systems are discussed. The conclusions are drawn in Section llVl In the Appendix[51 we give the algebraic 
cluster form factor expression for the four-nucleon projectile. 



II. CLUSTER FORM FACTOR AND SPECTROSCOPIC FACTOR CALCULATION 



In this section we derive expressions for the channel cluster form factors for a composite system of A nucleons, a 
projectile of a nucleons and a target of A — a nucleons. All the nuclei are assumed to be described by eigenstates 
of the NCSM effective Hamiltonians expanded in the HO basis with identical HO frequency and the same (for the 
eigenstates of the same parity) or differing by one unit of the HO excitation (for the eigenstates of opposite parity) 
definitions of the model space. We limit ourselves to a < 4 projectiles. In such a case, the projectiles can be efficiently 
described by a Jacobi-coordinate HO wave functions. The target and the composite system is assumed to be described 
by Slater determinant single-particle HO basis wave functions which is in general more efficient for A > 4. In this 
section we present results for a = 1,2,3. The cluster overlap integral for the a = 4 projectile is given in Appendix 
IbI The NCSM effective interaction theory is not repeated in this paper. It can be found in Ref. [3j for the case of 
two-nucleon interactions and in Ref. p4| for the case of two- plus three-nuclcon interactions. 



A. Coordinate and HO wave function transformations 

We follow the notation of Ref. [2f| . We consider nucleons with the mass m neglecting the difference between the 
proton and the neutron mass. For the purpose of the present paper we use the following set of Jacobi coordinates: 
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Here, £0 is proportional to the center of mass of the A-nucleon system: R = y^Co- On the other hand, £ p is 
proportional to the relative position of the p + 1st nucleon and the center of mass of the p nucleons. The f]A-a 
coordinate is proportional to the relative position between the center of masses of the two interacting clusters, i.e. 
the A — a nucleon target and the a-nucleon projectile. The 1? coordinates appear only for a > 1. Let us rewrite the 
Eq. and JTJi) as 
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where R,£ M a = [n + r 2 + ■ ■ ■ + f A - a ] and R a CM = yjl [r A - a +x + . . . + r A ] Following, e.g. Ref. 0, the HO 

wave functions depending on the coordinates J3J) transform as 

Y,(LMlm\Qq)ip NLM (RcM a ) l Pni m (RcM) = 

Mm 

{n'l'N'L'Q\NLnlQ)^(l'7n'L'M'\Qq) Vnn , ml {^ a ) VN , L , M ,{( ) , (3) 

n'l'm' N' L' M' 

where (n'l' N' L'Q\NLnlQ) » is the general HO bracket for two particles with mass ratio -^z^- 



B. Composite and asymptotic wave functions and the channel cluster form factor 

We consider the ^4-nucleon composite state eigenfunction 

(£ . . . ff A - a ■ ■ ■ $A-i<J X • ■ • OATl • . • T A \A\JMTM T ) (4) 

with the a and r the spin and isospin coordinates, respectively. The J and T is the total angular momentum and 
the total isospin, respectively, and M, Mt their third components. The A stands for the additional quantum numbers 
needed to characterize the eigenstate. The $ coordinates appear only for a > 1. 

Projectile-target wave function with the radial wave function describing the relative motion of the two nuclei 
replaced by the Dirac delta function can be written as 

ft t I * a Q i^(A-a,a)JMTM T x v 

(fl • • • U-a-lVA-a VA-a V A -a+X ■ ■ ■ "&A-X^X ■ ■ ■ &ATX ■ ■ ■ TA\^ aIlTl Jl 2 Tr,sl > "VA-al 

= J2(hMxI 2 M 2 \sm s )(sm s lmi\JM)(TxM Tl T 2 M T2 \TM T ) ^ (?M ~ a Z^A y tf A _ a ) 

^ ' VA-aVA-a 

x (£l • • ■ £,A-a-X<JX ■ ■ ■ ^A-aTx ■ ■ ■ TA-a\A - aaIxMiTiM Tl ) 

X 0A-a+X ■ ■ ■ $A-XVA-a+X ■ ■ ■ @ A T A—a+X ■ ■ ■ T A \af3I 2 M 2 T 2 M T2 ) , (5) 

where (£1 . . . ^A-a-x^T\A — aaliMxTxM^} and (^A-a+i ■ ■ ■ $ ^A-xa^aPIzMzT-iMT^) are the target and the projectile 
eigenstates, respectively. The s is the channel spin and the / is the channel relative orbital angular momentum. For 
our convenience, we also define a projectile-target wave functions with a HO radial wave function describing relative 
motion of the two nuclei, i.e. 

(Sl---SA-a-X VA-a ^A-a+l---VA-X(Jx---(TATx---TAW aIiTi J l2T2 . sl \ nl) 

= ^(/ 1 Af 1 / 2 M 2 | S m s )( S m s /TO ; |JM)(r 1 M Tl T 2 Af T2 |TM T )i? n; (^_ Q )y ;m! (r7 A _ Q ) 

x (£l • • ■ ^A-a-X^X • ■ ■ & A—aTX ■ ■ ■ T A -a\A ~ aaIxMxTxM Tl ) 

X i'&A-a+X ■ ■ ■ ^A-XVA-a+X ■ ■ ■ ^A^A-a+X ■ ■ ■ T A I a f3I 2 M% T 2 Mf 2 ) • (6) 



The R n i (r) in Eq. © is the radial HO wave function with the oscillator length parameter b = bo = y , where m 
is the nucleon mass. Due to our use of the coordinate transformations Q the oscillator length parameter is the same 
for all coordinates, i.e. &o- ln Eqs. ijoKil) . the coordinates $ appear only for a > 1. 
The channel cluster form factor can then be defined as 

ui X - J a T ahTl ,a0i 2 T 2 ;si(VA-a) = (AXJT\A^;% T T2;sl ;S VA _ a ) = ]T R nl (r, A - a) (AX JT\A^~^f^ sl ; nl) 



{A~a)\a\ 



with A the antisymmetrizer. As stated above, we assume identical HO frequency for all eigenstates and identical (or 
differing by a single HO excitation in the case of opposite parity states) definitions of the model space. 
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The spectroscopic factor is obtained by integrating the square of the cluster form factor. In particular, we have 

^A X - J aaI 1 T u a0I 2 T2;sl = Jdr]A-aVA-a\ u A- J aaI 1 T 1 ,af3I 2 T 2 \sl( r lA-a)\' 2 

= ^El(^l<w 2;s; ;«OI 2 . (8) 



As in this paper all the eigenstates are assumed to be expanded in a large but finite HO basis, we can set the 
integration limit to infinity in Eq. JHJ). 

It turns out that obtaining the eigenstates using the Jacobi coordinates becomes increasingly difficult with the 
number of nucleons A mostly due to the complicated antisymmetrization. As stated in the Introduction, for A > 4 
it is by far more efficient to use the SD basis. Consequently, it is desirable to to express the overlap Q using the 
eigenstates obtained in the SD basis. 

The relationship between the Jacobi coordinate and the SD eigenstates is 

(r*i . . .r A v\ ■ ■ .<taTi ■ ■ .T A \A\JMTM T ) S n = (£i . . .ff A - a ■ .-^a-i^i • • - oat\ ■ ■ .T A \A\JMTM T )tp 00 a(£ ) , (9) 

for the composite system and similarly for the A — a nucleon target. The subscript SD refers to the fact that this 
state was obtained in the Slater determinant basis, i.e. by using a shell model code, and, consequently, contains the 
spurious CM component. 

To arrive at the desired expression, we investigate an analogues overlap to using the SD eigenstates. We consider 
corresponding SD eigenstates to iQJ) and ©, i.e., 

sMAXJT\A^;^% T T2 . sl -,nl) S u , (10) 

where all the composite and the target eigenstate Jacobi coordinates are replaced by the Cartesian coordinates. The 
projectile eigenstate is kept unchanged with the Jacobi coordinates. Further, the rf A - a is replaced by -R^m- Explicitly, 
we have for the SD analog of the state JHJ 

(n . . . f A -a RcM $A-a+l ■ ■ ■ $A-lVl ■ ■ ■ O A T X ■ ■ ■ T A | ®tf~£*pi*£™ T > n 0sD 

= ^(/ 1 M 1 / 2 M 2 | S m s )( S m s ZmHJM)(T 1 A/ Tl T 2 M 

x (fi . . . r A ^ a ai . . . o a ~qT\ ■ ■ ■ T A - a \ A - aa/iMiTiM Tl ) SD 

X (#A-a+l • • ■ 3 A -l<T A -a+l ■ ■ ■ VATA-a+1 ■ ■ ■ T A \a(3I 2 M 2 T 2 M T2 ) . (11) 

We now proceed in two steps. First, using the relation @ for both the composite and the target eigenstate and the 
transformation @ we obtain 

S u{AXJT\A^;^ T2 . sl -,nl} so = (nl00l\00nll)^ (AXJTlA^^^nl) , (12) 

with a general HO bracket due to the CM motion, which value is simply given by 

2n + l 

(nl00l\00nll)_^_ = " ■ ( 13 ) 

The relation iff^l has been derived in the past, see e.g. Refs. 0], 0, [23] ■ Second, we relate the overlap flTTf to 
a linear combination of matrix elements of a creation operators between the target and the composite eigenstates 
s-D(AXJT\al lili j i . . . dn a i a j a \A — aaIiTi)su- The subscripts ni/iji refer to the single-particle state quantum numbers 
n i('i5)ji TO il m ti e t c - Such matrix elements are easily calculated by shell model codes. To obtain the channel cluster 
form factor we use the second equality in Eq. . 

C. Single- nucleon projectile 

In the case of a single-nucleon projectile, the asymptotic state (JHJ simplifies as no $ coordinates are present. 
The projectile wave function has just spin and isospin components with I 2 = \ and T 2 — i, respectively. It is 
straightforward to calculate the overlap of the states l@J and (JSJ) ■ The result is given by 

(AXJT\ A * ( ^" M) 1 T. ^6^) =VAj2Rni(VA-i)§j(-l) h+J+j \ | ") 

a'l-il, 22' — L u J J 

x {A\JT\{N A _ l i A _ 1 I 1 T l -nlj\)JT){N A _ 1 i A _ x I x T 1 \A-laI l T x ) , (14) 
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with s — \/2s + I. The composit e eig enstate is expanded in a basis with lower degree of antisymmetry using the 
coefficients of fractional parentage [2q 

((N A -iiA-ihTi',nlj^)JT\AXJT) = Y,N l ( N A-iiA-ihT 1 ;nlj\\NiJT){NiJT\A\JT) , (15) 

with N = iVX— i + 2n + I the total number of HO excitations for the A nucleons and i, «a-i the additional quantum 
numbers that characterize the A- and A — 1-nucleon antisymmetrized basis states, respectively. 

To obtain the overlap integral matrix element starting from SD composite and target eigenstates, we make use of 
HI 21 with a = 1 and perform the above discussed second step. That is quite straightforward for the a — 1 case and 
we easily arrive at the final expression: 

(a\jt\ a ^.i^--' - WioU A B£"'- iitw { 'i j ; 

x S D<AAJT|||a], Ij .|||A-laIiT 1 )s D . (16) 
Using Eq. JEJ, we obtain the familiar CM correction factor (^j)^ 0, IE 13, E3, E3, E3, 

D. Two-nucleon projectile 

For a > 1 projectiles we only present the overlap matrix elements for the composite and target wave functions 
obtained in the SD basis. For a — 2, which includes the deuteron projectile, the derivation is slightly more complicated 
due to additional re-couplings and explicit presence of the a — 2 relative coordinate wave function expanded in the 
HO basis. The final expression reads: 

(AXJT\ A ^lf^ sl -,8 VA _ 2 ) -E^'M-^ ^oOiloLfl , E^^la = 2f5I 2 T 2 ) 




x 55 f 5 ~iJ J 2 ( T\h+J+l+h+T 2 \ h h s 1 J / L afc Z 2 
x SS2 / 2 j a j b / ab i ab (-l) ( , J 7 a6 / \ S2 h hb 

x {n a l a n b l b L ab \nln 2 l 2 L ab ) 1 S B(AXJT\\\(al alaJa al blbn )^ T ^\\\A - 2aI 1 T 1 ) SB , (17) 

with the antisymmetry condition for the two-nucleon channels (— l) ; 2+ s 2+T 2 _ _j two-nucleon projectile wave 

function (n 2 l 2 s 2 I 2 T 2 \a = 2(3I 2 T 2 ) is expanded in the HO basis depending on $a-i- The spin and isospin components 
of the wave function depend on the spin and isopin coordinates cja-i&a, and ta_i,ta, respectively. For the deuteron 
projectile, I 2 = 1, T 2 = 0, s 2 = 1 and l 2 = or 2. Here, in addition to the HO bracket (|13[) due to the CM correction, 
one more HO bracket appears that corresponds to particles with mass ratio 1. This is due to the transformation of the 
HO wave functions (p n i m (R'^)ip n2 i 2m2 ('d a-i) to the single-particle HO wave functions (p na i a m a (rA)Vn b i b m b (rA-i)- 

E. Three-nucleon projectile 

For a = 3, which includes the triton or 3 Hc projectile, the derivation is still more complicated due to additional 
re-couplings and explicit presence of the a — 3 relative coordinate wave function expanded in the HO basis. The final 
expression reads: 

(AXJT\ A #itaf T2;s p^_ 3 ) = E^'fo-3) ^00i|00ntt)_,_ jf 

n a-3 

x ^{n 2 l 2 s 2 ht-i-M 2 £ 2 J 2 \)I 2 T 2 \a = 3/3J 2 T 2 ) 

X sIs 2 j2hj2jajbjcIab^ 2 Ll b (-l) 



2 T 2 , {\h-I+J+lc+l+J2+ I ab \ h k S 



I J I 



l a h L ab I f I A L 2 j 2 
1 1 s 2 W C 2 l c I ab I. 
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L 2 L ab l 2 

S2 J2 /ah J I 2 I " j I 7 1 • t 

( ja J6 J a6 J ^ J 2 2 Jc 

x (n a l a n b l b L ab \N 2 L 2 n 2 l 2 L ab )i(n c l c N 2 L 2 X\nlJ\f 2 C 2 \) 



x SD <AAJT|||((at a t J^)al, . )^ \\\ A - 3aI lTl ) su , (18) 
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The three-nucleon eigenstates are expanded in a basis with lower degree of antisymmetry using the coefficients of 
fractional parentage 25] 

{{n2hs2ht2;^2^J2\)hT 2 \a = 3f3I 2 T 2 ) = J2 m (n 2 l 2 s 2j2 t 2 ; N 2 C 2 J 2 \\\NiI 2 T 2 ){NiI 2 T 2 \a = 3f3I 2 T 2 ) , (19) 

with N = 2Af 2 + C 2 + 2n 2 + l 2 the total number of HO excitations for the three nucleons and i the additional 
quantum number that characterizes the three-nucleon antisymmetrized basis states. The 12 — j symbol of the first 
kind [23 appearing in Eq. 118(1 is defined in Appendix 1X1 For the triton or 3 He projectile, I 2 = \,T 2 — \ and 
(— l) J 2+ £ 2 _ 1 j n (|XH|). in addition to the HO bracket (|13fl due to the CM correction, two general HO brackets 
appear that correspond to particles with mass ratios 1 and i. These are due to the sequence of two transforma- 
tions of the HO wave functions ^ n im{RcM) ( Pn2i2rn 2 ('&A-i) ( PN'2C2M2{'^A- 2 ) to the single-particle HO wave functions 

(fA)yn b l b rn b {rA-l)<PnJ cmc {rA- 2 )- 



III. APPLICATIONS 



In this section, we present results of cluster form factor and/or spectroscopic factor calculations for ( 5 He| 4 Hc+n). 
( 5 He| 3 H+d), ( 6 Li| 4 He+d), ( 6 Be| 3 He+ 3 He), ( 7 Li| 4 He+ 3 H>, ( 7 Li| 6 Li+n), ( 8 Be| 6 Li+d), ( 8 Be| 7 Li+p), ( 9 Li| 8 Li+n) and 
( 13 C| 12 C+n). All calculations are done using the approach described in Sect. [H] The composite A-nucleon system 
and the target A — a nucleon system are described by the NCSM wave functions obtained in the m-scheme Slater- 
determinant basis shell model calculation. In particular, we use the many-fermion dynamics (MFD) shell model 
code and a specialized transition density code that calculates the (a) . . . a*) matrix elements employing the wave 
functions obtained by the MFD (or the Antoine |3|) code. The projectile a-nucleon NCSM wave functions for a = 3 are 
obtained in the Jacobi-coordinate HO basis using the code MANYEFF (2^. For the a = 2, i.e. the deuteron projectile, 
the relative-coordinate wave function is obtained using the standard NCSM two-body effective interaction code, see 
e.g. Ref. Q. As a technical point, we note that in the case of a = 2 there is no CM HO binding potential contrary 
the usual NCSM two-body effective interaction calculation. Consequently, the overlap of the full-space a — 2 wave 
function with the model s pac e P might not be large. This then could lead to numerical difficulties when applying the 
Lee-Suzuki procedure p9Tl30|| to obtain the model space a = 2 wave functions. To address this issue, higher precision 
than the double precision had to be used in the relevant computer code. 

We performed several calculations to test correctness of the formulas presented in Sect. |H]as well as their computer 
coding. First, we cross checked that the Eq. I|14f) and Eq. (|16f) give the same result for the ( 5 He| 4 He+n) system. In 
the former case, we employed the Jacobi-coordinate MANYEFF code while in the latter we used the SD basis MFD 
shell model code together with the transition density code. Obviously, the same effective Hamiltonian was used in 
both calculations. To test the Eqs. (|17I18|I we switched the role of the projectile and the target. For example, for 
the ( 3 H|d+n) system, we can apply Eq. (|16fl with the deuteron as the target and the neutron as the projectile or we 
can apply Eq. I|17|l with the neutron as the target and the deuteron as the projectile. In the latter case, the neutron 
target state is described as a single neutron N = HO state. Similarly, the Eq. (|18fl can be tested considering the 
( 4 He| 3 He+n) system described with the neutron as the projectile in Eq. I|16fl or 3 He as the projectile in Eq. JTSJl. 
Finally, we also performed a test for the ( 5 He| 3 H+d) system by switching the target and projectile and using the Eqs. 
(|17|l and l|18l) . respectively. We note that all these tests are non-trivial as the projectile and the target are described 
using different coordinates and the respective wave functions are obtained by different computer codes. 

In this Section, all the calculated channel cluster form factors are presented as a function of the separa- 
tion r between the CM of the projectile and the CM of the target. In particular, r — y jj^^ VA-a with 
ffA-a defined in Eq. ©)• Consequently, we have to use the reduced mass /i = ( A ^ a ) a m m ^he defini- 
tion of the HO length parameter, b = = \J (A-a)6p ° - presented channel cluster form factors are 
then related to those defined in Eq. Q by ^^^^(r) = £» i?„^r, 6)(AAJT|^$^ a j/^ 2;si ; nl) = 

X)„(^5r^) 3 ^ 2 -^ >ni ( 7 W- a ' bo)(A\JT\AQ^hTl a pi^T2-si-> "0- The spectroscopic factor (JSJ) is not affected by the choice 
of the coordinate: SfL™^^^ — J driA-aVA-a\ u A XJ a^i 1 T 1 ,a(3i2T2-,si('>lA-a)\ 2 - J ' drr 2 \ui x JJ ahTi af3l2T2 . sl {r)\ 2 = 
J2 n \{A^JT\AQ^j~^ a pj 2 T Tr s i', nl}\ 2 . The channel cluster form factors presented in this section are obtained from Eqs. 
(|16I17I18|I with the radial HO wave function R n i(rjA- a ,bo) replaced by R n i(r, b). 
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FIG. 1: (Color online) Overlap integral of the 5 He § ground state with the 4 He+n as a function of separation between the 
4 He and the neutron. Dependence on the basis size for JV max = 4, 6, 8, 10, 12 is presented. The CD-Bonn 2000 NN potential 
and the HO frequency of hQ, — 16 MeV were used. The s and I are the channel spin and the relative angular momentum, 
respectively. 



A. ( 5 He| 4 He+n) 



In Figs. njandGl we present our ( 5 He| 4 Hc+n) results for the 5 He | ground state resonance. The dependence of the 
channel cluster form factor on the basis size of the NCSM calculation is shown in Fig. ^for basis sizes from N meix = 4 
(4hil) to N max — 12 (12fiil). Here, N mSLX is the maximal number of the HO excitations above the unperturbed ground 
state. The CD-Bonn 2000 NN potential 3JJ and the HO frequency of Ml = 16 MeV were used in the calculations. 
Clearly, with increasing 7V max , the changes between successive curves become smaller, a sign of convergence. 

When calculating overlaps involving 4 He and a p-shell nucleus it is not obvious which HO frequency is the optimal 
one due to differences in the radii of the participating nuclei. In general, the HO frequency in the NCSM is typically 
fixed so that the binding energy has the least dependence on it. However, while for 4 He the fastest NCSM convergence 
and the least dependence on the HO frequency is obtained with a higher HO frequency, e.g. hfl > 20 MeV, for the 
p-shell nuclei the optimal frequency lies typically in the range of hfl = 10 — 15 MeV. In Fig. we present the channel 
cluster form factor dependence on the HO frequency using a wide range of frequencies: Ml = 13 — 19 MeV. It is 
satisfying that the sensitivity of the cluster form factor to the choice of the HO frequency is rather small. 

The same as in Figs. ^ and |5] is presented in Figs. |21 and 0] for the excited \ 5 He resonance. This resonance 
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FIG. 2: (Color online) Overlap integral of the He 4 ground state with the He+n as a function of separation between the 
4 He and the neutron. Dependence on the HO frequency for TlQ = 13, 16, 19 MeV is presented. The CD-Bonn 2000 NN potential 
and the basis size of iV max = 12 were used. 

is broader in experiment |3^|. Our calculation shows a more extended overlap integral for the i state compared 

to the | state. At the same time, both the basis size and the HO frequency dependencies are more pronounced. 

Nevertheless, even for this broad state the conclusions reached for the | state apply, which is an encouraging result. 

It can be seen in both Figs. ^ and [3] how with the increasing N max the overlap extends at large r. However, due to 
the finitcncss of our basis, the overlap integral approaches zero with increasing r even for states that correspond to 
physical resonances. 



B. ( 5 He|t+d) 



In Fig. \5\ we present the NCSM calculated channel cluster form factors for the ( 5 He|t+d) J w = | + resonance. The 
NCSM calculations were performed using the CD-Bonn 2000 NN potential in the 11M7 model space for 5 He. There 



are three possible channels corresponding to the combinations of the channel spin s = | and s 



and the relative 



orbital momenta I = and 1 — 2. Clearly, the overlap integral is by far the largest for the s-wave channel s = | 
/ = 0. The d-wave channels give small overlap integrals with the s = | I = 2 channel greater than the s = | / = 2 
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FIG. 3: (Color online) Overlap integral of the 5 He | first excited state with the 4 He+n as a function of separation between 
the 4 He and the neutron. Dependence on the basis size for 7V max = 4, 6, 8, fO, 12 is presented. The CD-Bonn 2000 NN potential 
and the HO frequency of HQ = 16 MeV were used. 



channel. Results corresponding to three different HO frequencies are shown in Fig. Despite the wide range of the 
HO frequencies, changes in the resulting overlap integrals are not significant. 

It is interesting to point out that the d+t | + resonance in 5 He appears as the second f + state in the NCSM 
calculations reported in this paper as well as in the previously published 5 He results (3^, EU- The appearance of 
low- lying positive parity states in the NCSM calculations of the 5 He spectrum was criticized in Ref. |3j|- We note 
that such states were also observed in phenomenological shell model calculations [3|||37|]. At the same time, some 
evidence for a low- lying ^ state was reported in the R-matrix analysis of Ref. (3^. No such state was, however, 

included in the recent evaluation |32j |. The low- lying | + state obtained in the present calculation has basically zero 
overlap with the d+t. It is quite possible that it corresponds to a non-resonant continuum state of a free neutron and 
the 4 He that appears as an excited state due to the fmiteness of the basis used in our investigation. The d+t resonant 
| + state is dominated by the s 3 p 2 configuration. 
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FIG. 4: (Color online) Overlap integral of the 5 He i first excited state with the 4 He+n as a function of separation between 
the 4 He and the neutron. Dependence on the HO frequency for hQ, = 13, 16, 19 MeV is presented. The CD-Bonn 2000 NN 
potential and the basis size of iV max = 12 were used. 

C. ( 6 Li| 4 He+d) 

Our ( 6 Li| 4 He+d) channel cluster form factors for the 6 Li J 71 " = 1 + , 3 + and 2 + are presented in Fig. |SJ Fig. [TJand Fig. 
|H1 respectively. The corresponding spectroscopic factors are then summarized in Table [I] In the NCSM calculations, 
we used the 6 Li wave functions obtained in Ref. [3,9] using the CD-Bonn NN potential |40j. In the three figures, the 
thick lines correspond to the 1QMI results and the thin lines to the 8fiJ7 6 Li results. We can see only small changes 
in the overlap integrals when the basis size is changed, in particular for the l + ground state and the 3 + excited 
state. It is interesting to note that the ground state is dominated by the s = 1, I — 4 Hc+d configuration while the 
excited 3 + 0, 2 + and 1^0 states are dominated by the s = 1, I = 2 4 Hc+d configuration. This is in agreement with 
the analysis of the 4 He+d elastic scattering experimental data [32| . 



D. fBe | 3 He+ 3 He) 



Just as a resonance plays a critical role in the rate of the d+t reaction, there is some speculation that the reaction 
3 He( 3 He,2p) 4 He, which is important to the Standard Solar Model (SSM), could be affected by a resonance in the 
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FIG. 5: (Color online) Overlap integral of the 5 He | + excited state with the 3 H+d as a function of separation between the 3 H 
and the deuteron. Dependence on the HO frequency for hQ = 13, 16, 19 MeV is presented. The CD-Bonn 2000 NN potential 
and the basis size of -/V max = 12 (for 3 H and d) N max = 11 (for 5 He) were used. 



composite 6 Be system |4l|. Although recent experiments at the LUNA underground facility [42j seem not to favor 
a narrow resonance, they do not definitively rule out its presence |4lj . To investigate a possibility of a resonance in 
the 3 He+ 3 He system, we performed NCSM calculations of the overlap integrals ( 6 Be| 3 He+ 3 He) for the lowest four 
+ l states obtained in the NCSM description of 6 Be. The calculations were performed using the CD-Bonn 2000 NN 
potential in the basis spaces up to lOTitl for 6 Be. The lowest two + l states are the p-shell dominated states while the 
third and the fourth + l state is a one-particle-one-hole and a two-particle- two- hole dominated state, respectively. 
In the N max — 10 (lQhfl) basis space and the hfl = 13 MeV calculation, their excitation energy is 12.5 MeV and 
13.5 MeV, respectively, not far from the 3 He+ 3 He threshold. However, the excitation energy of these states is not 
yet converged in the present calculation and it is expected to further decrease with the basis size enlargement. Our 
channel cluster form factor results obtained in the 10fi.fi space are shown in Fig. |5| The dependence on the HO 
frequency is presented for all four states. A large overlap integral is found for the ground state and also for the 
one-particle-one-hole dominated 0^1 state. On the other hand, the overlap integral for the 0^1 state is negligible 
and the one for the 0^~1 is quite small. It is interesting to note a stronger HO frequency dependence of the overlap 
integrals for the 2hil dominated states compared to the p-shell states. This is another manifestation of a slower 
convergence of these states in the NCSM. The significant overlap integral of the 0^1 6 Be state suggest that this state 
might contribute as a resonance in the 3 He+ 3 He reaction. However, our prediction of its excitation energy is not 
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FIG. 6: (Color online) Overlap integral of the 6 Li 1^"0 ground state and the 1^0 first excited with the 4 He+d as a function of 
separation between the 4 He and the deuteron. Dependence on the basis size for Nmax = 8, 10 is presented. The CD-Bonn NN 
potential and the HO frequency of TiQ. = 13 MeV were used. 



certain. Based on our NCSM results up to N max — 10 we expect this state to converge below the 3 He+ 3 He threshold 
of 11.49 MeV. 



E. ( 7 Li| 4 He+t) 

Our results for the ( 7 Li| 4 Hc+t) channel cluster form factors are shown in Fig. ^3 while the corresponding spec- 
troscopic factors are summarized in Table ITT1 Apart from the large overlap integrals and spectroscopic factors for the 
bound | x and 5 states we find these quantities large also for the first excited | and the first excited | state. 
Both these states appear as resonances in the 4 He+t cross section [sjj. The present results can be compared to the 
three-nucleon transfer calculations of Ref. 0] obtained using the phenomeological Cohen-Kurath interaction [43| . 

The agreement for the lowest four states is quite good. For the second excited | state, however, our spectroscopic 
factor is significantly smaller than the one obtained in Ref. [l9j . 
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FIG. 7: (Color online) Overlap integral of the 6 Li 3 + first excited state with the 4 He+d as a function of separation between 



the He and the deuteron. Dependence on the basis size for N u 
HO frequency of TiQ, = 13 MeV were used. 



10 is presented. The CD-Bonn NN potential and the 



F. ( 7 Li| 6 Li+n) 

The other system involving 7 Li as the composite nucleus that we investigated is 6 Li+n. Our calculated overlap 
integrals are summarized in Fig. 1111 The corresponding spectroscopic factors are given in Table GUI As in the 
( 7 Li| 4 He+t) case, we observe large overlap integrals and spectroscopic factors for the two bound states | a and ^ . 
Contrary to the ( 7 Li| 4 He+t) case, however, we find a large overlap integral and the spectroscopic factor for the | 
state. The lowest | x and jL states have negligible overlap integrals for the 6 Li+n system. The large overlap integral 
and the spectroscopic factor for the | 2 state is consistent with the observed resonance in the 6 Li+n cross section. In 
addition to the | state, we also find large overlap integrals for the higher lying | and | states. In Fig. we 
display the basis size dependence of the and the | states for the iV max = 4, 6 and 8 calculations. The results for 
the resonant | state are fairly robust. The spectroscopic factor of the | state show a stronger basis size dependence. 

Our N max = 8 spectroscopic factor for the jL state is 0.016 for s = |, I = 1 channel, a significantly smaller value 
than that obtained using the Cohen-Kurath Op-shell phenomenological interaction [l5|]. Our NCSM result is, however, 
consistent with the spectroscopic factor obtained using the variational Monte Carlo wave functions |44j . 
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FIG. 8: (Color online) Overlap integral of the 6 Li 2 + excited state with the 4 He+d as a function of separation between the 
4 He and the deuteron. Dependence on the basis size for N m ax = 8, 10 is presented. The CD-Bonn NN potential and the HO 
frequency of fiQ. — 13 MeV were used. 



G. ( 8 Be| 6 Li+d) 



We also investigate systems with 8 Be as the composite nucleus. The 6 Li+d reactions in particular are of some 
interest in controlled thermonuclear research and their cross sections have been measured j4£j. At the same time, 
they are a part of a reaction network with 8 Be as the composite nucleus, which is being analyzed by the R-matrix 
method |46|. Our calculated spectroscopic factors for the 6 Li+d channels are presented in Table llVl The J n = 2 + 
channel cluster form factors are then shown in Fig. ^| We label the J 71 ' 8 Be states as they are obtained in the 
current 67ifi calculation. The description of the excitation spectra of 8 Be in the NCSM is generally very good ^tJ- I n 
Ref. 47] in addition to the p-shell states, the slowly converging intruder + 0, 2+0 and 4+0 states were found. Such 
states have complicated structure with wave functions dominated by higher than OhQ components. Existence of such 
states was controversial 0, 0, 0, E3, 13 • However, in the latest evaluation a broad intruder 2 + state 
is included at about 9 MeV excitation energy. Such a state is required by the R-matrix fits of nuclear reactions that 
involve 8 Be as the composite system. While in Ref. the intruder states were investigated in the basis spaces up 
to 10?iO, in this paper we use the 6h£l wave functions to calculate the channel cluster form factors. In this space 
the intruder states appear at a higher excitation energy and their importance is likely suppressed because of that. In 
the present calculations, the + intruder state is the state 0^" and the 2 + intruder state is the state 2^". Even in 
the current 6fi.fi basis space, these states have significant overlaps with the 6 Li+d system. We note that the 6 Li+d 
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TABLE I: Spectroscopic factors for the ( 6 Li| 4 He+d) corresponding to the 6 Li ground and excited states and the 4 He ground 
state. The CD-Bonn NN potential, the basis size of iV max = 10 for 6 Li and iV max = 12 (for 4 He and d) and the HO frequency 
of hQ = 13 MeV were used. The s and / are the channel spin and the relative angular momentum, respectively. 
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FIG. 9: (Color online) Overlap integral of the 6 Be 0+0 states with the 3 He+ 3 He as a function of separation between the 3 He 
nuclei. Dependence on the HO frequency for hi} = 13, 16 MeV is presented. The CD-Bonn 2000 NN potential and the basis 
size of iV m ax = 10 (for 6 Be), iV max = 12 (for 3 He) were used. The 0+1 and 0+1 are p-shell states. The 0+1 and 0+1 are a 
one-particle-one-hole dominated and a two-particle-two-hole dominated state, respectively. 
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FIG. 10: (Color online) Overlap integral of the 7 Li low-lying J = | , | , | , | states with the 4 He+ 3 H as a function of 
separation between the 4 He and the triton. The CD-Bonn 2000 NN potential, the basis size of iV m ax = 8 (for 7 Li), AT max = 10 
(for 4 He and 3 H) and the HO frequency of HQ = 13 MeV were used. 



has a rather high threshold of 22.28 MeV. Interestingly, we obtain a dominant overlap integral in the s = 2, I = 
channel for the 8 Bc 2 + excited state number seven, which is a p-shell O7if2-dominated state in our calculation with 
the excitation energy of 22.54 MeV. 

We note a technical issue affecting the ( 8 Be| 6 Li+d) overlap integral calculations using Eq. I|17(l . When an m-scheme 
calculation is performed with a fixed M, it is in general necessary to use M > for J > channels in order to generate 
all needed reduced matrix elements of the operators. 



H. ( 8 Be| 7 Li+p> 

Our spectroscopic factors for the ( 8 Be| 7 Li+p) channels are presented in Table [Vj The overlap integrals for the 
J n = 1 + channels are shown in Fig. It should be noted that the expressions for cluster form factors and 

spectroscopic factors presented in Sect. ^ employ isospin formalism. In order to distinguish, e.g. a proton or a 
neutron projectile, the Eq. H16|l must be multiplied by the isospin Clebsch-Gordan coefficient (TxMt^MtJTMt) 
with Mt 2 = +5 ( — h) for proton (neutron) j5||. This coefficient is 1 for all the overlaps studied in this paper except 
the ( 8 Be| 7 Li+p) overlap for which it is equal to — 4=, -4= for T — 0, 1 8 Be states, respectively. We find large overlap 
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TABLE II: Spectroscopic factors for the ( 7 Li| 4 He+ 3 H) corresponding to the 7 Li ground and excited states and the 4 He ground 
state. The CD-Bonn 2000 NN potential, the basis size of iV max = 8 (for 7 Li), 7V max = 10 (for 4 He and 3 H) and the HO frequency 
of HQ = 13 MeV were used. The s and I are the channel spin and the relative angular momentum, respectively. 
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TABLE III: Spectroscopic factors for the ( 7 Li| 6 Li+n) corresponding to the 7 Li ground and excited states and the 6 Li ground 
state. The CD-Bonn 2000 NN potential, the basis size of iV max = 8 and the HO frequency of TiQ. = 13 MeV were used. The s 
and I are the channel spin and the relative angular momentum, respectively. 

integrals for the if and the lj states as well as large spectroscopic factors for the 3^~ state consistent with the 
resonances in the 7 Li+p cross section |3 E^- We also note very large spectroscopic factors and cluster overlap 
integrals for the lj state which is the second T = 1 1 + state in our calculation appearing at the excitation energy of 
20.37 MeV. Such a state is not included in the current evaluations 0,15^. It is, however, needed in the R- matrix 
analysis [46j. Comparing our spectroscopic factors to the Cohen-Kurath calculations [15|. we have a reasonable 
agreement for the lowest states (after correcting for the above discussed factor of 2 due to the isospin Clebsch-Gordan 
coefficient), with increasing differences for the higher lying states. Our results are also influenced by the isospin mixing 
of the 2J, 2^ and if, 1% states. 

I. ( 9 Li| 8 Li+n> 

The experimental information on 9 Li is rather limited l57| . New experiments are under way or planned, 
however, to explore this nucleus. One of such experiments is the inverse-kinematic d( 8 Li,p) scattering 58] . It is 
therefore useful to perform theoretical calculations of the ( 9 Li| 8 Li+n) spectroscopic factors. Our results obtained for 
both the negative and the positive parity states of 9 Li are summarized in Table IVII In the present calculation we 
employed the AV8'+TM'(99) two- plus three-nucleon interaction in the NCSM calculations performed in an approach 
described in Ref. j^. The Tucson-Melbourne (TM) three-nucleon interaction was introduced in Ref. E3 with the 



18 



0.3 



0.25 - 



1 0.2 



* 0.15 

<u 

-i— 1 



0.1 



0.05 



/ 

I 

I ; 
I : 
I : 
I : 

I 

I •' 
I : 

I:// 

.7 

1:1 



7' 



i 

i y 



V 

•/ 

oh 



< 7 Li I 6 Li+n > 

\ /=l/2"3/2"5/2"7/2" 



\ 

\ 

. \ 
•. \ 

•. \ 

*• \ 

\ 

•• \ 



\ 



\ 



1/2/= 


1 


3/2 j 


3/2 /= 


1 


3/2" 2 


1/2 /= 


1 


\I2\ 


3/2 /= 


1 


l/2" 2 


3/2 /= 


1 


5/2" x 


3/2 /= 


1 


5/2 2 


1/2 /= 


3 


112 \ 



_1_ 







4 

r [fm] 



FIG. 11: (Color online) Overlap integral of the 7 Li low-lying J = | , | , | ,2 states with the 6 Li+n as a function of 
separation between the 6 Li and the neutron. The CD-Bonn 2000 NN potential, the basis size of iV max = 8 and the HO 
frequency of fiQ. = 13 MeV were used. 



particular version we are using, the TM'(99), described in Ref. [6(]]. The Argonne V8' NN potential is a slightly 
simplified version of the high-quality AV18 interaction [lj. Our current calculations are limited to the 4M7 and 5hQ 
basis spaces for the negative and the positive parity states, respectively, due to the complexity of the calculation with 
a genuine three-nucleon interaction. In general, it is accepted that a three-nucleon interaction is needed in addition to 
the high-quality NN potentials to explain the few-nucleon system binding energies and to improve description of some 
three-nucleon scattering observables. Recently, it has been shown that the three-nucleon interaction is also needed 
for a correct description of low- lying excitation spectra of p-shell nuclei Q, |24| . 

In Fig. 1151 we show the channel cluster form factor for the lowest two | states. The | x state with a large 
spectroscopic factor is a candidate for the 4.296 MeV 9 Li state lying just above the 4.063 MeV 6 Li+n threshold. 



J. ( 13 C| 12 C+n) 

Apart from an increase of the binding energy, the genuine thrcc-nucleon interaction also causes an increase of 
the spin-orbit splitting. This is demonstrated not only in different level spacing and sometimes in a different level 
ordering in calculations with the three-nucleon interaction, e.g. in 10 B, but also in the spectroscopic factors and the 
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FIG. 12: (Color online) Overlap integral of the 7 Li low-lying | states with the 6 Li+n as a function of separation between 
the 6 Li and the neutron. Dependence on the basis size for iV max = 4, 6, 8 is presented. The CD-Bonn 2000 NN potential and 
the HO frequency of hQ. — 13 MeV were used. 



overlap integrals as we show in this subsection for ( 13 C| 12 C+n) system. Using the wave functions obtained in Ref. 

|24j . we compare in Fig. Eland Table IVIll the h and | channel cluster form factors and the spectroscopic factors, 
respectively, obtained in calculations with and without the TM'(99) three-nucleon interaction. The AV8' NN potential 
is used for the two-nucleon interaction. We can see that the | channel cluster form factor and the spectroscopic 

factor increase when the three-nucleon interaction is included, while at the same time the § factors decrease. This 
can be understood as an increase of the spin-orbit splitting of the Op | -Op i levels due to the three-nucleon interaction 
which results in a purer (Opi) 8 12 C ground state and a purer (Opa) 8 , (Opi) 1 13 C ground state. 

Comparing to the phenomenological Cohen-Kurath spectroscopic factors |15| , a better agreement is achieved in the 
more realistic calculation with the three-nucleon interaction. 



IV. CONCLUSIONS 



We derived expressions for calculations of channel cluster form factors and spectroscopic factors from the ab initio 
no-core shell model wave functions. We considered the most practical case, with the composite system and the target 
nucleus described in the Slater determinant harmonic oscillator basis while the projectile eigenstate expanded in 
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TABLE IV: Spectroscopic factors for the ( s Be| 6 Li+d) corresponding to the 8 Be ground and excited states and the 6 Li ground 
state. The CD-Bonn 2000 NN potential, the basis size of iV max = 6 and the HO frequency of Ti£l = 13 MeV were used. The s 
and I are the channel spin and the relative angular momentum, respectively. The 0+0 and 2+0 are intruder states. All other 
are p-shell states. 

the Jacobi coordinate HO basis. The spurious center of mass components present in the SD bases were removed 
exactly. The calculated cluster form factors are then translationally invariant. The algebraic expressions for the 
channel cluster form factors were derived for up to four-nucleon projectiles. We numerically tested these expressions 
for systems consisting of up to a three-nucleon projectile. Several numerical tests were performed that involved 
interchanging the role of the target and the projectile as well as performance of two independent calculations one of 
which employed only the Jacobi-coordinate wave functions for all nuclei involved while the other used the SD basis 
wave functions for the composite system and the target. Identical results were obtained in both cases. 

As examples of application, we presented results for ( 5 He| 4 He+n), ( 5 He| 3 H+d), ( 6 Li| 4 He+d), ( 6 Be| 3 He+ 3 He), 
( 7 Li| 4 He+ 3 H), ( 7 Li| 6 Li+n), ( 8 Be| 6 Li+d), ( 8 Be| 7 Li+p), ( 9 Li| 8 Li+n) and ( 13 C| 12 C+n) systems, with all the nuclei 
described by multi-?if2 NCSM wave functions. The calculations involve no fitting. Apart from the basis size, the only 
parameter appearing in the NCSM is the HO frequency, which is typically fixed so that the binding energy is the least 
dependent on the HO frequency. In the current application, this is hard to achieve in the cases that involve both the 
0s- and the Op-shell nuclei, as we require the same HO frequency for all nuclei. Therefore, we studied the dependence 
on the basis size and on the HO frequency in most investigated cases. It is very encouraging that our results ere 
rather stable and robust. Additionally, we found a qualitative agreement with experiment for, e.g., ( 5 He| 4 He+n), 
( 5 He| 3 H+d), ( 7 Li| 4 He+ 3 H), ( 7 Li| 6 Li+n), ( 8 Be| 7 Li+p), in the sense that large channel cluster form factors correspond 
to resonances in cross sections. This confirms that the multi-frfi NCSM wave functions provide a realistic description 
of light nuclei, in particular for the low- lying p-shell states. 

As a further development, apart from performing calculations for systems with a four-nucleon projectile, our goal is 
to utilize the channel cluster form factors as a first step to describe low energy reactions on light nuclei. The presently 
calculated channel cluster form factors were obtained using the model space wave functions. As the next step, we 
need to take into account the influence of the complementary space and calculate effective, or renormalized, channel 
cluster form factors. It is expected that this will improve the cluster form factors at intermediate distances and make 
them more suitable for matching to the correct asymptotic cluster wave functions. Hopefully, it will be possible to 
develop a microscopic nuclear reaction approach similar to the RGM |J| starting, however, from realistic ah initio 
wave functions. 
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FIG. 13: (Color online) Overlap integral of the 8 Be 2 + states with the 6 Li+d as a function of separation between the 6 Li and 
the deuteron. The CD-Bonn 2000 NN potential, the basis size of -/V ma x = 6 and the HO frequency of TiQ. = 13 MeV were used. 
The 2^~0 is an intruder state. All other are p-shell states. 
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APPENDIX A: 12 - j SYMBOL DEFINITION 

The 12 — j symbol of the first kind is defined by 
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FIG. 14: (Color online) Overlap integral of the 8 Be 1 + states with the 7 Li+p as a function of separation between the 7 Li and 
the proton. The CD-Bonn 2000 NN potential, the basis size of iV mm = 6 and the HO frequency of HQ — 13 MeV were used. 
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TABLE V: Spectroscopic factors for the ( 8 Be| 7 Li+p) corresponding to the 8 Be ground and excited states and the 7 Li ground 
state. The CD-Bonn 2000 NN potential, the basis size of iV m ax = 6 and the HO frequency of Ti£l = 13 MeV were used. The s 
and / are the channel spin and the relative angular momentum, respectively. The 0+0 and 2+0 are intruder states. All other 
are p-shell states. 



APPENDIX B: FOUR-NUCLEON PROJECTILE 

The channel cluster form factor expression for the case of a four-nucleon projectile is 
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TABLE VI: Spectroscopic factors for the ( 9 Li| 8 Li+n) corresponding to the 9 Li ground and excited states and the 8 Li 2 + l 
ground state. The AV8'+TM'(99) two- plus three-body interaction, the basis size of Af max = 4,5 and the HO frequency of 
HQ = 14 MeV were used. The s and I are the channel spin and the relative angular momentum, respectively. Only channels 
with S > 0.001 are shown. 
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TABLE VII: Spectroscopic factors for the ( 13 C| 12 C+n) corresponding to the 13 C ground and excited states and the 12 C ground 
state. Results obtained using the AV8'+TM'(99) two- plus three-body interaction and the AV8' NN interaction are compared. 
The basis size of iV max = 4 and the HO frequency of hQ, = 15 MeV were used. The s and I wee the channel spin and the relative 
angular momentum, respectively. All the presented 13 C states are p-shell states. 



Similarly as the three- nucleon eigenstates in Eqs. (|18I19|I . the four-nucleon eig enstates are expanded in a basis with 
lower degree of antisymmetry using the coefficients of fractional parentage |25j | 

{{{{n2l2S 2 j2t2;N2C 2 J 2 \)J i T z )N^J s \)hT 2 \a = A(3I 2 T 2 ) = Y^{ n ^ s ^^^ C ^k\\ N ^ J ^) 

x (N x i x J 3 T 3 ;Af 3 C 3 J 3 l)\\NiI 2 T 2 ) 
x (NihT 2 \a = ApI 2 T 2 ) , (B2) 

with N x = 2M2 + C 2 + 2n 2 + l 2 and N = N x + 2J\f 3 + C 3 the total number of HO excitations for the three and 
four nucleons and i x ,i the additional quantum numbers that characterize the three- and four-nucleon antisym- 
metrized basis states, respectively. In the case of the 4 He projectile, I 2 = T 2 = and (— i) l 2+^2+c 3 _ \ j n 
Eq. l)Bl|l . in addition to the HO bracket (jlHfl due to the CM correction, three general HO brackets appear that 
correspond to particles with mass ratios 1, ^ and ^. These are due to the sequence of three transformations of the HO 

wave functions ip n i m (R^)ip n2 i 2m2 (&A-i) W 2 £ 2 m 2 (^-2) W 3 c 3 m 3 to the single-particle HO wave functions 

^nj a m a (rA)^n b i b rn b (^A^i)'Pnj cma {^A- 2 )^n d i d rn d ( ff A-3)- With the help of Eq. (jBlJ, one can study the alpha-cluster 
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FIG. 15: (Color online) Overlap integral of the 9 Li two lowest | states with the 8 Li+n as a function of separation between 
the 8 Li and the neutron. The AV8'+TM'(99) two- plus three-nucleon interaction, the basis size of iV max = 4 and the HO 
frequency of hfl — 14 MeV were used. 
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FIG. 16: (Color online) Overlap integral of the 13 C \ ground state and the | first excited state with the 12 C+n as a function 
of separation between the 12 C and the neutron. Results obtained using the AV8' NN potential and the AV8'+TM'(99) two- 
plus three-nucleon interaction are compared. The basis size of iV max = 4 and the HO frequency of TiQ. — 15 MeV were used. 
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structure of the p-shell nuclei states. In the past, this has been typically investigated using the cluster models, see 
e.g. Ref. 0. 
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